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Abstract 

The set of common roots of a finite set I (it is an ideal) of homogeneous polyno- 
mials is known as projective algebraic set V . In this article I show how to dualize 
such projective algebraic sets V by elimination of variables from a system of poly- 
nomials with the Grobner bases method. A dualization algorithm is implemented in 
Singular sec GPS09 . Some examples are given. The main diagram shows the 
relationship between the ideal I, its radical yl and their dual ideals. 

1 Introduction 

1 read about the duality principle in the book |Gie| and saw some examples of du- 
alization in the books |Gie| . |Kom| for quadrics and in the more recent introductory 
book on plane curves |Fis) for plane curves. This last book gives some particular 
examples how the dualization is carried out but no general method. Some authors 
mention that variables have to be eliminated from a system. For plane curves the 
system is derived nicely in |Bri] . For one hypersurface I recently found [PW] p. 104f. 
But for intersections of hypersurfaces the only example that I found was the inter- 
section of two hypersurfaces in |Kom| as a general example, without being specific. 
So I derived in this article the system for the intersection case. 

In what follows I recommend to read |GP| for the theoretical background about 
projective space, homogeneous polynomial, ideal, projective variety, etc., which is 
not covered in the present work. For an introduction to Grobner bases see |Cox| or 
[Fro] . I used the methods derived in this article to dualize some examples and also 
checked with the examples given in |MP| . 

2 Motivation with plane curves 

In what follows here, we assume that the denominators do not vanish. Think of the 
inversion radius r as having value i = (Other values are also permitted, e.g. 

1.) We consider different representations of plane curves. 
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2.1 Parametric 

We consider a parametrically given plane curve c(t) — (x(t), y{t)Y . Then we can 
define the pedal curve of c(t) with respect to the origin as 

_ y(t)x'(t)-x(t) y'(t) ( -y'(t) 
m)) ~ x'itf + yHtf \ At) 

(the pedal is the locus of the feet of perpendiculars from the origin to the tangents 
of the curve c(t)). We can also define what it means to invert c(t) with respect to 
the circle of radius r around the origin: 



X(t) 

x(t) 2 + y(t) 2 V. y(t) 



By composing the two maps given above we get the dual of c(t) as the inverse of the 
pedal: 



d{c(t)) = i(p(c(t))) = 



-y'(t) 



y(t)x'(t)-x(t)y'(t) V. x'(t) 
This is best explained by a commutative diagram: 



c(t) 



Y 




L(p( C (t)))^^p(c(t)) 



2.2 Complex 

Now we do the same for a curve z(t) in the complex plane. The pedal is 

2z'(t) 

The inverse is 

i(z(t)) = 4= 

z(t) 

and the dual is (again by composition) 

2r 2 z'(t) 



d(z{t)) = t(p(z(t))) = 



z(t) z'(t) - z'(t)z(t) 
A similar commutative diagram as in the parametric case holds. 

2.3 Implicit 

For implicitly given curves f(x, y) = we cannot give explicit formulas for the pedal 
curve but we give a method for computing it. We need the gradient Vf(x,y) — 

(jg£(x,y), ^(x,y)^j . Assume p — (x,y) is a point of / and P — (X, Y) is a point of 
the pedal of /. Then, by the definition of the pedal, following must hold: 
1. p = (x, y) is a point of /: f(x, y) = 0. 
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2. P lies on the tangent at / through p: (P — p)* V/(x, y) = 0. 

3. P is orthogonal to tangent at / through p: P* y), — §^(£, J/)^ =0 

By eliminating (a;, j/) from these three equations we get an equation in (X, Y) which 
is the pedal curve. For convenience we substitute (X , Y) h-> (x,y). In what follows, 
we will see how the elimination can be done with Grobner bases. 

The inverse of f(x, y) = is / f x 2 + x y 2 1 ^s+yi ) = 0- The dual of / is the compo- 
sition of inversion and pedal as constructed above. 



3 Theory 

3.1 Case of one homogeneous polynomial 

First we consider the following projective algebraic set 

V(p) = {x G K n+1 I p(x) = 0} 

with K an algebraic closed field, x = (xo, Xi, ■ ■ . , x n ) a point of K™ +1 and p a 
homogeneous polynomial from K[a;o, Xi, . . . , x„]. V(p) consists of all roots x = 
(xo, Xi, . . . , x n ) of p and is a hypersurface in the projective space P n . 

Let u = (uo, Ui, ... , u n ) be a normal vector to V(p) in a regular point x G V(p). 
On the other side we know that the gradient 

'*)=(s w * H S H 

is normal to V(p) in x. Therefore u and Vp(x) are linearly dependent. This can be 
written as u = AVp(x) with a factor A. We can form the following system 

p(x) = 

(1) 

u - AVp(x) =0 

Definition 3.1.1. We define the setV{p) = {u G K n+1 j p(x) = 0, u-AVp(x) = 0} 
0} partial solutions to the system ^ to be the dual ofV(p). 

Note that we are not interested in a complete solution of but only in the 
partial solutions which I call here the u-part of the solution. The u-part of the 
solution of this system is the result of applying the Gaufi map to p(x) = 0, where the 
Gaufi map (see 0, p. 103) is given - in Chow coordinates (see [L]) - by 

7 : x h-> u = AVp(x) 

Now we want to construct an equivalent system to |l]) but simpler in structure, 
describing the same dual algebraic set V*(p). 

Theorem 3.1.2. There exists a system B of polynomials in K[wo, wi, • ■ • , u n ] with 
the same solution set V*(p) as the system £7]). 

Proof. We start with the system viewed as system of polynomials 
qj G K[a;o,xi, • ■ ■ ,x„, X,u , «i, . . . , u n \ 
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and eliminate the first n + 2 variables. The system is: 

qi{x , xi,... ,x n , A, u 0l U!, . . . , u n ) = p(x , . . . ,x n ) = 
?2 (%0, Xi, ■ ■ ■ , X n , A, Wo, Ml 



dp 



L, • ■ • , u, nj 



Uo — A- (Xo, Xn) = 

OX 



dp 



qi+2(xo, Xi, . . . ,X n , A, Wo, Ml, ■ ■ ■ ,U n ) — Ui — A- — (x , ...,£„) = 

OXi 



q n +2(xo,Xl, . . . ,X n ,X,Uo,Ul, . . . ,U n ) = U n — A- (x , ■ ■ ■ ,X n ) = 

OX n 

Let G be a Grobner basis for the ideal S — (qi, . . . , qj, . . . , qn+3) with respect to an 
elimination ordering, where e.g. x > A > u. By the Elimination Theorem of |Cox| 
this basis G eliminates x and A and B = G n K[«o, Ui, . . . , u n ] is a Grobner basis of 
the elimination ideal E — S n K[uo, Ui, . . . , u n ]. We get E = (B). By construction, 
this basis B is a system of polynomials in K[uo, ui, . . . , u n ] having the same partial 
solution set V*(p) as the original system □ 

Because of this property we define E = (B) to be the dual ideal of the initial ideal 
(p). With the canonical isomorphism Ui Xi we can map the ideal E G K[u] to an 
ideal D £ K[x] which leads to the following commutative diagram: 



(p) € K[x 



D G K x 




E e Kful 



Example 1. As an example we dualize the quadric (n = 3): 



l 2 . 01 bi — 1 2 . 1 a , 1 2 , 
-box + ZbiXoXi Xx + —x 2 + —x 3 = 



a-i 



a-z 



a.', 



The system is here: 



—box + 2b\x§x\ 



6l — I2.I2.I2 n 

X\ H x 2 H 2:3 = 

01 a 2 a 3 

26o2;oAi — 2biXi\\ + no = 

26j - 2 



-261X0A 



-x\ Ai + «i = 



«2 



-X2A1 +M2 = 



a; 3 Ai + u 3 = 

03 

By eliminating Xi and A /rom these equations we get one equation in m, which is the 
dual quadric 



(61 — 1)wq + 2ai6iwo«i + ai£>o«i + (aioafei — a2&obi + a 2 6o)'«2 + (010361 — 03^061 + 0360)113 = 
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Example 2. As another example we dualize a quadric over Q[w,x,y, z]: 



(w 2 -x 2 ) 

and get 




How about dualizing this result? This simple example shows that we must be able to 
dualize ideals given by more than one polynomial. This is done in the next section. 



3.2 Case of a system of homogeneous polynomials 

The same method can be extended and applied to finite sets of polynomials which 
represent geometrically their intersection as hypersurfaces. We consider the projec- 
tive algebraic set 



v ■=v(i) = v( Pl , p a ,...,p fn ) = {xe 



r n+l 



|Pi(x) 



: P m (x) = 0} 



with K an algebraically closed field, x = (xo, xi, . . . , x n ) a point of K. n+ and / the 
ideal generated by m > 1 homogeneous polynomials pi,p2, ■ ■ ■ ,Pm from K[x]. 

The dual projective algebraic set V* of V is also a projective algebraic set which 
is the zero set of an ideal E generated by polynomials from K[tto, tti, . . . , u„]. We get 
this elimination ideal E by eliminating (using a suitable Grobner basis) the m + n+1 
variables Xi and Xj from the m + n + 1 equations: 

Pi(x , ...,x n ) = 



Mo 



p m (x , ...,x n ) = 
^\j-£^{x ,. . . ,x n ) =0 



2 = 1 



m - ^ Xj^-(x , ...,x n ) = 



(2) 



j'=i 



-5> 



dx 



' 3 f Xo, ■ ■ ■ ,X n ) = 



2 = 1 



dx 



The Xi are point and the m hyperplane coordinates but their roles can be inter- 
changed. 

With the Jacobi matrix of p(x) = (pi,P2, • • • ,Pm)*(x): 



/ dpi dpi 

dxQ dx\ 



J = 



dpi \ 

dx n ' 
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and setting A = (Ai, A2, . . . , A m )* we can write the system ^ vectorially as 

p(x) = 

771 

u V" AjVpj(x) = u - J*(x)A = 

3=1 

For m = 1 this is the same as system We redefine the Gaufi map: 

7 : x H-> u = J (x)A 

Example 3. We dualize now the following i 

2 2 

Z2 - Z3 5 
£0 — £2 

TTie system |S|) is here (written as ideal) 



(3) 



/ 22 



\ 



/ 



Xo — SE2, 

— A 2 + uo, 

Mi, 

— 22: 2 Ai + A 2 + «2 

\ 2x 3 Ai + u 3 
A sorted Grobner basis with the elimination property is 

( u i> \ 

Uq + 2tioll2 +u\—u%, 

A2 — Uo, 
2x 3 Ai + u 3 , 
X2W3 + a^tto + X3U2, 
X2U0 + X2U2 + X3U3, 

2X2 Al — Uo — U2, 

2 2 
2:2-2:3, 

\ 2:0-2:2 / 

and we see that the first two elements generate the elimination ideal E for this example 

Ui, 

Uq + 2U()U2 + u% — u\ 

3.3 Main diagram 

We have, by denoting with D(yT) the dual of the radical ideal of I: 




D(y/I) 



y 1 

V(D(VI)) v(^Jd(V7)) V(D(I)) V(^D(T)) 
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(The bent arrow in the middle of the diagram needs a proof.) As an example ex- 
plaining this diagram consider: 



I = 

Then we have 



x -\-y — z 



D(VJ) = \Jd(VJ) = ( x - y ) 
D(I), 



x-y, 
y 2 + 2yz + z 2 

y + z, 

x-y 



VW)-- 

By adding one polynom we get another example: 

( ^ 
I = \ x + y- z, 



D(I) = 



D{VT) = sfoiyT) = {x-y) 
xy - y 2 + xz - yz, \ _ f (x- y)(y + z) 



x 2 -2xy + y 2 J \ (x-y) 2 
^W) = { x-y ) 



4 Examples 

4.1 Steiner's Roman surface 

I wrote a Singular procedure dual using the fast elimination provided by Singular 
with the combination of hilb and eliminate for computing the dual ideal D of a given 
ideal /. The procedure dual takes as argument a finitely generated homogeneous ideal 
I — (pi, . . . ,p m ), where all generating polynomials pt are homogeneous and elements 
of the basering. 

A first example Singular input file for testing the functionality is given in the 
appendix. This file dualizes Steiner's Roman surface (see MP for other nice dual- 
ization examples): 

Example 4. 1. First we set the basering to Q{xq,x 1 ,X2,Xz\. 

2. As an example we want to dualize the ideal 

I=( x\x\ - xqXxx 2 xo, + x\x\ + x\x\ ) 

3. We call procedure dual on this ideal and the first thing which it does is to adjoin 
auxiliary variables to our ring. The new ring is: 

Q[x , xi , x 2 , xz , Ai , u , ui , u 2 , ti 3 ] 
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4- The next step is constructing the following ideal generated by the system W of 
equations, which in this case is 

/ 2 2 ,22,22 \ 

/ x x x 2 - x x 1 x 2 x 3 + x 1 x 3 + x 2 x 3 , \ 

XiX 2 X S \i + Uq, 
— 2XlX 2 \l + X X 2 X3\! — 2x 1 x%\ 1 + Ul, 

—2x\x 2 X\ + a;oa;iX3Ai — 2x 2 x^Xi + u 2 , 
\ x x 1 x 2 \ 1 — 2x(x' i X 1 — 2x\x- A Ai + u 3 / 

5. The Grobner basis for S w.r.t. the elimination order is (we only show some 
elements at the beginning and the end here because it is bigger than this page) 

I 4u% — U u\ — W m1 + U!U 2 U[i — U0U3, ^ 

+ X3A1M0M2 — xlXl^U^i + UoU!U 2 — 2uqU 3 , 
_ -ixlM'ul — xlXxu's + U U!U 2 — 2u%U3, 

XqX 2 Xz\i — 4xiX 2 x^Xi — 4x 2 xsXi — 2xoa;ia;3Ai + xoui + 2x 2 u 3 
\ a^^i^Ai — ixiXsXi — 4x\x 2 X3Xi — 2xqx 2 X3Xi + xqu 2 + 2x\U3 J 

6. The elimination ideal E consists here only of the first element of G 

E — ( 4mq — uqUi — uou 2 + U1U2U3 — U0U3 ) 

This is the dual ideal of I. We can interpret the x as point coordinates and the 
u as hyperplane coordinates. But we want to be able to pass this ideal E to dual 
and dualize it too! For this to work, we have to make one final step in dual and 
this is to map the xs to the us and viceversa, giving as result the dual ideal of 
I in point coordinates, which is 

D = ( 4xq — X x'f — X X 2 + X\X 2 X3 — xqx\ ) 

7. When calling dual on D we get I as dual ideal of D. We won't show here the 
steps, but the reader can generate them using the example input file from the 
appendix. 

The input file to Singular for this example generates the figure [7] by using 
surf ex. lib which in turn uses the program surf. 

4.2 Parametrized quadric 

We present now another example, the forth and back equalization of a somewhat 
special quadric. 

Example 5. 1 . Our ring is at the beginning Q(ai, a 2 , 03, bo, bi)[xo, Xi, x%, Xz\. 

2. The quadric is given by the zero set of the following ideal (this zero set is a 
projective algebraic set) 

Q = ( -bozo + 2&iZoxi - bj ^xl + ±xl 

3. We adjoin auxiliary variables and the new ring is 

<Ql(ai,a 2 ,a3,ba,bi)[xo,x 1 ,X2,X3,X 1 ,u ,u 1 ,u 2 ,U3] 
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Figure 1: Steiner surface and its dual surface 



4- Here we construct the system |3] of equations 



I -bo^o + 2bix x 1 



S = 



b -^xi + ±xl \ 



V 



2b x >^i — 2b 1 x 1 \ 1 + u , 

-2&lXoAl + ^^iClAl +M1, 
- — X2X1 + U2, 
U3 



5. The Grobner basis of S is 
( 



G = 



"3. 



(&i — 1)mq + 2ai&iii iti + ai&o - "? + (a 1 a 2 bl — a 2 b b 1 + a 2 b )u 2 , 
2x 2 \i — a 2 u 2 , 
{aia 2 b\ — a 2 bobi + a 2 bo)xiu 2 — aibix 2 uo — 01602:2^1, 
(2ai6j — 26 ( )&i + 2b )x 1 Xi — aibxuo — aib ^i, 
a 2 boXou 2 — a 2 b\X\u 2 + x 2 uo, 
aiboxoui — (b± — 1)xi_uq — 2aibix±ui ~ a\b\X 2 u 2 , 
Xotto + x x u x + x 2 u 2 , 
26 s Ai — 2&1X1A1 + u , 
aia 2 boXo — 2aia 2 biXoXi + (a 2 bi — a 2 )x\ — a\x\ 



6. The elimination ideal ( the dual of Q) is 



E : 



(61 — 1)u'q + 2ai&iuoiii + aibouf + (aia 2 b\ — a 2 bobi + a 2 bo)u 2 , 
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7. We substitute it h-> x in E and get 

,.2 



D 



(61 — l)x + 2aib 1 x x 1 + a 1 b x 1 + (aia 2 b 1 — a 2 b b 1 + a 2 bo)x 2 , 

X3 



This ideal D can be now again dualized which constitutes our next example. 
Example 6. 1. After adjoining the auxiliary variables, our ring is 

Q(ai,a 2 ,a i ,bo,b 1 )[x ,x 1 ,x 2 ,X3,\ 1 ,\ 2 ,u ,ui,u 2 ,U3] 

Notice here that - since our ideal D has two generators - we have to adjoin two 
As (the procedures do this automatically) . 

2. The system Q - in fact an ideal too - is here 

( (61 — \)x\ + 2a\b 1 x x 1 + ai&oXi + (a^bl — a 2 b bi + a 2 b )x%, \ 

X3, 

— (2&i — 2)a; Ai — 2ai6i^iAi + u , 
—2aibiXa\i — 2aiboXi\± + mi, 
— {2a\a 2 b\ — 2a 2 bobi + 2a 2 bo)x 2 \i + u 2 , 
— \ 2 + U-J, 



s 2 = 



V 

3. The correspon 
( 



Grobner basis is 



G 2 = 



a\a 2 bau — 2aia 2 biUoUi + (a 2 bi — 02)1*1 — a\u 2 , 
M — «3, 

X3, 

(2aia 2 bl — 2a 2 bob\ + 2a 2 bo)x 2 \i — u 2 , 
a\x\u 2 — a\a 2 bix 2 uo + (a 2 bi — a 2 )x 2 u\, 
(2a\b\ — 2aibobi + 2aibo)xi\i — a\b\uo + (&i — l)ui, 
xou 2 + a 2 box 2 uo — a 2 b\x 2 u\, 
(b\ — l)xoui — aiboxiuo + 2aibiX\ui + aibix 2 u 2 , 
X0U0 + X1U1 + x 2 u 2 , 
(2&i — 2)a; Ai + 2aibi:riAi — u , 
\ (61 — 1)xq + 2aibiXoXi + axboxj + (aia 2 &? — a 2 b bi + a 2 b )x 2 ! ) 



4- The new elimination ideal is 

( a\a 2 boUQ — 2aia 2 &iu iti + (a 2 b 1 — a 2 )u\ — aiw| ) 
5. After mapping the xs and us we finally get Q again 

( aia 2 b()xl — 2a\a 2 bixoXi + (a 2 bi — a 2 )x\ — a\x\ ) 

One might wonder why this looks different than the initial equation for Q, but 
if you multiply the original equation with the non-zero factor —a 1 a 2 7^ - an 
operation which does not change the zero set, you see that you get the same 
quadric. 
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4.3 8-shaped space curve 

Example 7. In this example we intersect a sphere and a cylinder to get an 8-shaped 
space curve, which is given by the ideal: 

x 2 +y 2 ~l, 
x 2 + y 2 + z 2 - 2x - I 

You can see these surfaces and their intersection in figure^ 

Now we dualize the ideal I and get the following ideal (after dehomogenizing) with 
one polynomial as generator: 

I 4x 6 + 12x 4 y 2 + 12x 2 y 4 + 4y 6 - 12x 4 z 2 - 24x 2 y 2 z 2 - 12y 4 z 2 - 15x 2 z 4 + ' 
D = V2y 2 z 4 - 4z 6 + 36s 3 z 2 + 36xy 2 z 2 + 18xz 4 - 8x 4 - 16x 2 y 2 - 8y 4 - 20x 2 z 2 ~ 
\ 20y 2 z 2 + z 4 - Axz 2 + 4x 2 + Ay 2 

D is the dual of the 8-shaped space curve. They are depicted in figure [5[ 




Figure 2: 8-shaped curve as intersection of cylinder and sphere 

The Singular code for this example is: 

// load libs 

LIB "duality. lib"; 

LIB "surf ex. lib" ; 

ring r=0 , (t ,x,y ,z) ,dp; // ring over Q 
short =0; // print polynomials with 

// two polynomials 
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Figure 3: 8-shaped curve and its dual surface 



poly cylinder = x~2+y~2-l; 

poly sphere = (x-1) ~2+y~2+z~2-2~2 ; 

// intersection ideal (a space curved shaped like an 8) 

ideal il = cylinder, sphere; // inhomogeneous for plotting 

ideal i2 = homog(il, t) ; // homogeneous for dualising 

ideal dl = dual(i2); // dual of intersection ideal (a surface) 

poly d2 = subst(dl[l], t, 1); // dehomogenize dual for plotting 

d2; // show result 

// plot everything 

// (in surfex you may want to set transparency options for some surfaces) 
plotRotatedList (list (cylinder , sphere, il, d2) , list (x,y ,z) ) ; 

4.4 Examples from the introductory book [FisJ 

As further examples that the procedure dual works, I dualize some planar algebraic 
curves from the book Fis without including the intermediary output from Singular. 

Example 8. The Neil parabola is given by 

( X 1 - X X 2 j 

and its dual is 

( 4x1 + 27x x 2 2 ) 
You can see an illustration in figure [^] 
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Figure 4: Neil parabola (red) and its dual 

Example 9. The Newton knot is given by 

( x xl + x( — x x% ) 

and its dual cardioid is 

( 4x xf - 4xf + 27xqx 2 - 36x X!xl + %x\x\ - 4xi ) 
You can see an illustration in figure [5| 




Figure 5: Newton knot (red) and its dual cardioid 
Example 10. The hypocycloid is given by 

(00 o o o \ 

X Xi — 2x XiX2 — 2xoX 1 X2 + x x 2 — 2xqXix 2 + x 1 x 2 ) 

and its dual is calculated as 

( Xo + 3xqXi + 3xoxl + x\ + 3xqX 2 — 2\xqx\x 2 + 3x\x 2 + 3xqx 2 + 3x\x% + x 2 ) 

You can see an illustration in figure [6| 

Example 11. The next example was inspired by the Klein quartic. I changed one of 
the ellipses to a hyperbola and the generated ideal is 

I „4 5 2 2,1 4 3_2 2 , 15 „2 2 14 12 \ 
\ x ~ 4 X X 1 ~r 4 X 1 ~ 4 X X 2 -f j^XiX 2 — jX 2 70 / 
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Figure 6: Hypocycloid (red) and its dual 



and its dual ideal is calculated by Singular as generated by the polynomial 

I 123908754 2 - 1202643754° 4 + 44299185044 - 8228080004a;? + \ 
82762848044 - 43182720044° + 918881284 2 + 21866254°4- 
148231125444 + 902043450444 - 1921126200444+ 
1725988320444 - 5607878404° 4 - 11645585044 + 713525750444- 
784988540444 - 703298400444 + 91453593644+ 
23214240044 - 539359800444 - 507564960444- 
59801472044 - 19768672044 + 58816800444 + 11916134444+ 
\ 8018304044° + 3272256044° - 127534084 2 / 

You can see an illustration in figure ^ 




Figure 7: A quartic (red) inspired by Klein's quartic and the dual curve 
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5 Appendix 

The code for the Steiner surface example is: 
///////////////////////////////////////////// 

// This procedure calculates the dual ideal of the homogeneous ideal id 
// The output is a homogeneous ideal in the same ring 
///////////////////////////////////////////// 
proc dual (ideal I) { 
def RO=basering; 

if (npars(R0)>0) {ERRORO'Use a base ring without parameters!");}-; 
if (ord_test(R0) != 1 ) {ERROR (" The base ring must have a global ordering!");}; 
if (homog(I) !=1) {ERRORO'The input ideal must be homogeneous!");}; 
// get some information about the base ring and the input ideal 
int n=nvars(R0); 
int m=ncols(I) ; 

// change variables and compute transposed Jacobi matrix of I 

def NR=changevar("x() " ,R0) ; 

setring NR; 

ideal I=fetch(R0, I) ; 

matrix J=transpose(jacob(I)) ; 
// adjoin auxiliary variables to the ring 

def El=extendr ing (m, "10 " , "dp" , 1 ,NR) ; 

def R=extendring(n ) "u()","dp" ) l,El) ; 

setring R; 

matrix J=fetch(NR, J); 
// set up system S 

ideal I=fetch(NR, I) ; 

matrix L=matrix( [u(l . .n)] )-J*matrix( [1 (1 . .m)] ) ; 
ideal S=I,L; 

// eliminate first m+n variables from S by Groebner bases method 
int j , k ; 
poly prod=l; 

f or (k=l ; k<=n ; k++) {pr od=pr od*x (k) ; } ; 

for(j=l;j<=m;j++) {pr od=pr od*l ( j ) ; } ; 

intvec v=hilb(std(S) , 1) ; 

ideal Il=eliminate(S,prod,v) ; 
// resubstitute variables, such that the output can be used again as input 

map f=R, (u(l. .n) ,1(1. .m) ,x(l. .n)) ; 

ideal I2=ideal(f (II) ) ; 
// restore initial ring and return the result 

setring RO; export RO; 

return(fetch(R, 12)); 

> 

///////////////////////////////////////////// 
LIB "surf ex. lib" ; // Load library 
ring Rl = , (x(0 . . 3) ) ,dp; // Define ring 
// First example: Steiner 's Roman surface 

ideal I = (x(l) *x(2) ) ~2+(x(l) *x(3) ) ~2+(x(2) *x(3) ) ~2-x(0) *x(l) *x(2) *x(3) ; 
ideal D = dual (I) ; 
ideal DD = dual(D) ; 
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I; D; DD; // show results 
// Plot with surf ex 
ring R2 = 0, (x,y,z) ,dp; 
map f=(Rl, l,x,y,z); 

plotRotatedList (list (f (I) , f(D)), list (x,y ,z) ) ; 
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